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1 Introduction 



Quantum fields with impurities (defects) are recently subject of intensive studies in the 
context of condensed matter physics [1]- [5], conformal field theory [6,7] and integrable 
systems [8]- [16]. Most of the papers concern the case of one point-like defect. For more 
realistic applications however, the generalization to n > 1 defects would be essential. In 
the present paper we face the multi-defect problem by extending the algebraic technique 
developed in [14]- [18] for n = 1. We will show below that this technique perfectly applies 
to the case n > 1 as well, the key point being the reformulation of the boundary value 
problem at hand in algebraic terms. 

An immediate application of this formalism is the derivation of the Casimir force 
acting on a test particle. Since this issue is of some experimental interest, we treat it in 
detail, deriving the general expression of force. We consider also a complex scalar field and 
determine the relative charge density. As an explicit example we focus on two delta-type 
defects. It appears that the intensity and direction of the force depend on the position. 
The force is discontinuous at the defects and actually diverges when one is approaching 
any of them from the left (right). A compensation occurs however in the symmetric limit, 
which is finite. 

The paper is organized as follows. Introducing the notation, in the next section we 
briefly summarize the algebraic approach to one defect. This approach is generalized in 
section [31 where the solution for two defects is established. The latter is extended to 
n > 2 defects in section HI In section [3] we establish the two-point correlation function 
both at zero and finite temperature. Using this results we compute the Casimir force 
(section l5.2p and the charge density of the complex scalar field (section 15.41) . In section 
[6] we consider in details two delta-defects and discover a kind of resonance effect in the 
behavior of the density between the two defects. Some remarks and our conclusions are 
collected in section [71 

2 One defect: reflection-transmission algebra 

It is useful to recall first the algebraic treatment [15] of one defect on the real line. We 
take a defect with arbitrary position y G M. It divides the real line in two domains 
-Di = {x G M : X < y} and D2 = {x E W : x > y}. For simplicity we focus on the real 
massive scalar field (^{t,x), satisfying 

[d^-dl + m']ip{t,x)=0, x^y. (2.1) 

It is well known [19] that all possible dissipationless interactions of '^{t, x) with the defect 
y are described by the boundary condition 

( ititfy) ) = ( c ) ( dM^-y) ) ' ^ ^ ' ^'-'^ 
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where 

ad-bc=l, a,...,deR. (2.3) 

Following [17,18], one can solve equations fl2.1|2.2p by introducing the so called reflection- 
transmission algebra with the following structure. With each domain Di one associates a 
"creation" and "annihilation" operator {a*{k), ai{k) : k G M} which satisfy 

hik) , aj{p)] = = [a* (A;) , a*{p)] , 

[ai{k) , a*{p)] = [6,j+Tij{k)]27T6{k-p) + Rij{k)27T6{k + p), (2.4) 

where 

represent the reflection and transmission coefficients from the impurity. In terms of the 
parameters {a, ...,d} one has 

p &A;^ + i(a-6;)/c + c a% ^ 2iA; 

^^^^^^ = 6F + i(a + ^)fc-c ^ ' ^''^^^ = bk^ + iia + d)k-c ^ ^'-^^ 

^^^^^^ - &F-i(a + d)A;-c' ^^^^^^ " - i(a + rf)A; - c ' ' ^^'^^ 

Let us stress that the reflection coefficients (but not the transmission ones) depend on the 
position y of the impurity: in the following, we will also use the 'bare' coefficients Rjj{k) 
which are just the part independent from the position y: 

R,,{k) = Rii{k) e-^'^y , R22{k) = R22{k) e-^'^y . (2.8) 

The associated defect scattering matrix 

S{k) = T{k) + R{k) (2.9) 

satisfies unitarity and Hermitian analyticity 

s{k)s{ky = I , s{ky = s{-k) . (2.10) 

Now, the solution of (12. ip can be written in the form [17] 

- P [a*(A;)e'"('^)*-*'=" + a,(fc)e-^"('^)*+''="l . (2.11) 

-oo 27r^2a;(fc) 

where u{k) = \/k^~+lv?. The boundary condition (12. 2p then implies that {a*(A;), ai{k)} 
and {02(^)1 o,2{k)} are not independent, but satisfy the constraints 
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a^ik) = ^(i?.,(fc)a,(-/.)+T,,(A;)a,(A;)), (2.12) 
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which are consistent with fl2.4p because of (12.1 01) . The algebra 2li defined by equations 
(12 .4^ 12.12^ I2.13P is known as refiection-transmission (RT) algebra [15]. Introducing the 
matrices 

equations (12.12^ 12.131) can be conveniently written in the form 

A{k) = C{k)A{-k), A{ky = A{-ky C{ky . (2.15) 

By construction the matrix C{k) satisfies 

C{k)C{-k) = 1 , C{ky = C{-ky, (2.16) 

where the apex t denotes transposition. 

We end this summary with two observations. In the domain of {a, ...,d} where bound 
states are absent [17], the field (fi{t,x) satisfies in addition to (12.11 12. 2p the equal-time 
canonical commutation relations 

[ip(t,xi) , ip{t,X2)] = , [{dtip){t,xi) , ip(t,X2)] = -iS{xi - X2) . (2.17) 

We will not consider in this paper the issue of bound states, referring for the details about 
their impact to [17]. 

We stress in conclusion that the above construction is purely algebraic and applies 
for any representation of the RT algebra 2ti. This fact allows to investigate the physical 
effects associated with nonequivalent representations and in particular, to introduce [17] 
the concept of temperature, which is relevant for the applications in section 5. 



3 The RT algebra of two defects 

We show now how the RT algebraic framework allows to compute in a very systematic 
and straightforward way the global refiection and transmission coefficients when there 
are two defects. We confirm the results of [20], where these coefficients where computed 
through the calculation of Green functions. 

We consider two defects located at yi < y2- At each of them 



a) \ I 0,a Oct 



a = l,2. (3.1) 



Let R^°'^{k) and T^°'^(k) be the refiection and transition matrices defined by (12. 7p in terms 
of the parameters {oq, do}. According to (12.140 we set 



k) R^Si-k) r ^ M r|l(-fc) Rfli-k) 
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which satisfy by construction 

C^"\k)C^"\-k) =1 , C^"\ky = C^°'\-ky , a = 1,2. (3.3) 

The two defects separate now the real hne in three domains 

Di = {x eR : X < yi} , D2 = {x e R : yi < x < 1/2} , D3 = {x e M : y2 < x} . 

(3.4) 

Following the previous section, we associate with each of them the operators {a*{k), ai{k)} 
and define the field if{t,x) by (12. lip , where now i = 1,2,3. Since ip{t,x) must satisfy 
(13. ip . we impose 

A^''\k) = C^^\k) A^^\-k) , A^^\k)^ = A^''\-ky C^^\k)^ , a = 1,2, (3.5) 
where 

The operators {a* (A;), aj(/i;)} describe the propagation of (/?(t, x) in Dj. Only{ai(^), ai(fc)} 
and {a'^{k), a^{k)} give origin to asymptotic states. Therefore it is natural to expect that 
{a2{k), a2{k)} can be expressed in terms of {a\{k), ai(fc)} and {03 (/c), a^{k)}. Indeed, 
from (13.51) it follows that 

Cfl{k)a2{~k) = a2{k) - Cfl{k)a,{k) , 
Cg(-A;)a2(-A;) = a2{k) - C'il{-k)ai{k) . (3.7) 

Eliminating a2{—k) from these two equations, one gets 

'cfKk) - ctK-k)] a2{k) = CfKk) c'iK-k) a,{k) - Cg(-fc) Cfl{k) a,{k) . (3.8) 

Remark that because of the factors e^''^^^ and e^''^^^ appearing in i?22(^) -^22 (^) 
spectively, the coefficient Cf}{k) — (—/;;) never vanishes for generic 2/1,2/2 and (13.81) 
uniquely determines 02 (A;) in terms of ai{k) and 03 (A;). One has 

41(fc)rff(-fe)ai(fc) - R^kk)TS{k)a3{k) 
R^iik) - R^iik) 

Analogously, a2{k) is given in terms of a*(/c) and a^ik), 

^ 41(-fe)7^l'2'(-fe)«l(fc) - 4'j(-fe)7j2'(fe)«3(fe) .3 
i?i(-fc)-4^1(-^0 
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We conclude therefore that in order to define the RT algebra for two impurities 2I2, it is 
enough to fix the commutation relations between {a* (A;), ai(/c)} and {a^{k), 03 (A;)}. For 
this purpose we have to determine the global refiection and transmission matrices 



Ti3(A:) 
%i{k) 



7^(A;) 



7^ll(A:) 
7^33(A;) 



(3.11) 



in terms of the local ones and r[°]. According to fETT^ . r{k) and n{k) can be read 
from the relations between ai{k) and 03 (/c), which can be derived as follows. Plugging 
the expression (13. 8p into (13.51) . we get 



i-cSik)cr:ik) 



[2], 



ai{-k) + 



C12 (^) ^ui^) 



a^ik). (3.12) 



By means of the identities 

^22 i~k) ^ — ^22 (k) + ^22 i~k) ^ ^12 (k) ^21 i~k) 

^12 (~^) C*22 (— A;) ^ = —C[}{k) ^ c|2(A;) , 
following from (13. 3p . equation (I3.12p can be rewritten in the form 



ai{k) 



Cl5(A;)Cfl(A:)Cfl(fc) 



l-Cg(fc)Ci?(fc) 
The same type of calculations provides 



^12 (^) ^12 (^) 
- C22 (A;) C[}{k) 



[2], 



ai{-k) + 



a3(A;). (3.13) 



a 



[2], 
22 > 



-k) 



C21 (~^) C'L' (~^) C'22(— A:) 



l-CSi-k)C?}{-k) 
Comparing (13.131 131141) with (I2.12p . one infers 

nnik) 
nn{k) 



asi-k) 



a 



mi 
21 1 



k)d^}i-k) 



l-CS{-k)C?}{-k) 



cSik) 



dSik)cSik)cSik) 



i-cg(^)cfM) 

TS{k)Ti'^{-k)R^kk)e'''y^^ 



7^n(A;) e'^''^?^^ 



2]/ 



7^33(A:) 
^33(A;) 



a 



22 ' 



-fc) 



a 



[2], 



21 ' 



^) ^12 (~^) C'22 (~^) 



i-cg(-fc)cltx-fc) 



[2], 



7^33(A:)e-2i'=^^ 



[2],.^ , Ti^ik)Ti^i-k)R^Sik)e-'-^'y^^ 



1-B!ikk) 4'^(-A:)e-2ifc!'i2 



ai{k) . 
(3.14) 

(3.15) 
(3.16) 



(3.17) 
(3.18) 
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1 - eg {k) Ci? {k) 1 - 41 i-k) R^i (k) e^i^s/- ' 

1 - (-^) CS i-k) 1 - (/c) 41 e-2i%- ' 

which express the global two-impurity reflection and transmission coefficients in terms 
of the local ones. As for the one-impurity case, we have introduced 'bare' reflection 
coefficients which, as the transmission coefficients, depend only on the distance between 
the two impurities, yi2 = yi — y2- 

Remark that when one of the two defect is 'trivial', i.e. Ru{k) = R22{k) = and 
Ti2{k) = T2i{k) = 1, the global coefficients TZ and T become identical to the ones of the 
remaining defect. 

The relative S'-matrix is 

s{k) = r{k) + n{k). (3.21) 

Using the expressions (^JTMM) and the properties of R^°'\k) and T^'^^k) one can verify 
that fl3.2ip satisfies both unitarity and Hermitian analyticity (12.101) . which are fundamen- 
tal for the construction of the RT algebra 2I2 below. 

Equations fl3.17tl3.20|) describe in a compact way all (possible infinite) processes of 
reflection and transmission from the two defects, relating an incoming to an outgoing 
wave. Accordingly, we impose on the 2l2-generators 

[ai{k) , aj{p)] = = [a* (A;) , a*{p)] , 

[ai{k) , a*{p)] = [6i,+%j{k)]27rS{k-p)+nij{k)27r6{k + p), (3.22) 

where i,j take only the values 1 and 3. The exchange relations involving {a2{k), 02 (/c)} 
are determined by 03.221) and the relations fl3.9l I3.10p . 

Summarizing, the RT algebra 2I2 for two defects has the following simple structure. 
The elements (A;), 02 (Ac)}, associated with the field evolution between the defects, are 
linear combinations of {a* (A;), ai{k)} and {a^{k), 03 (A;)}, which generate the asymptotic 
states and satisfy quadratic exchange relations (13.221) . The latter have the same form as 
in the one-impurity case (12. 4p . but involve the global two- impurity reflection and trans- 
mission matrices Tl{k) and T{k), which are uniquely determined by the single defects 
via (13. 15113. 2DI) . We will show in the next section that this result has a straightforward 
generalization to an arbitrary number > 2 of defects. 
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4 General case of several defects 

We consider here n defects on M with coordinates yi < ■ ■ ■ < i/n- They define the n + 1 
domains 

Di = {x eR : X < yi} , Di = {x eR : yi-i < x < yi} , Ai+i = {x e R : yn < x} . 

(4.1) 

The relative RT algebra 2t„ is therefore generated by {a*{k), ai{k) : i = l,...,n + 1}. 
Extrapolating from the previous section, we expect that all {a*{k), ai{k) : i = 2,...,n} 
can be expressed in terms of {al{k), ai{k)} and {a^_^_i{k), a„+i(/c)}, which generate in 
turn the asymptotic states. Let us now verify this scenario and determine the algebra 21^- 
Setting 

(A;j= , a-l,...,n, (4.2) 

\ -^(0+1)0 1 -"'(a+l)(a+l)l ) 

we infer from the boundary condition (13. ip that (13.51 13.91) hold now for any a = 1, ...,n. 

Indeed, applying recursively the technics used in the case of two defects, one deduces 
that 

= [^if^^^^ , j,A; = l,...,n + l, (4.3) 

obeys the relations 

J^if^(k) = C[^-'](A;)AW1(-A;), (4.4) 
C^-A^ky = C^^-^\-ky and C^^-^\k) C^^-^\-k) =1, l<j <l<n + l ,(4.5) 

where C^^-^\k) is some matrix built on the C^°'^{k) matrices, a = j, . . . ,1. Starting from 

a,ik) = Ct'\k)a,{k) + d^,-^\k)a,{-k), (4.6) 
a,(fc) = ct''^'\k)aj{-k) + C^,-''^'\k)an+,{k), (4.7) 

one gets 



C^^r^'^{k)-C^^,-^\k)) a,{k) = C'-^\k) C^,-^\k) a,{k)-C'-=\k) Cf^^^k) a^^,{k). 

(4.8) 

Therefore, the expressions found in the case of two defects also apply here, with the rules 

CW(A;) C^'-^\k) ; C^^\k) ^ C^'-''+^\k) , 
ai{k) -> ai{k) ; a2{k) aj{k) ; a^{k) a„+i(A;) . 

Hence, it remains to compute the matrices C^^'''^\k) from the matrices C^°'\k). We give 
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here a recursion formula for them. We consider the relations 

AW(A;) = CW(A;)ylW(-fc) , (4.9) 

= CW(-A;)* and CW(A;) CW(-fc) = I , (4.10) 

^[2,n+i]^^) = C[2-"+il(A;)y4P'"+i](-A;), (4.11) 

which amounts to take j = 2 in the previous calculations. This leads to 

^[i,n+i](^) = c1^-"+'1(A;)A[^'"+i1(-A;), (4.12) 

or in components 

ai{k) = c!l-"+'l(^)ai(-fc)+Ca-""'"(^)«n+i(-^), (4.13) 

a„+i(fc) = c[l-"+'l(A^)ai(-fc)+c£-"'""(^)«n+i(-fc), (4.14) 



where 



" i-cS{k)crr^^k) ' ^'-''^ 



_ Cfi-"^^^(A;)cSkfc) .^x 



" i-cg(i^)cfr"+^i(-fc) • ^'-''^ 

Hence, an iterative use of these formulae, starting from the 'local' coefficients C["](/c), 
allows to compute all the C^^---^\k). 

Alternative recursion formulas can be obtained particularizing the last defect instead 
of the first one. One gets 



cl-\k) ct,^'\k) 

-^^(fc)a;-"'(fc) 

c£-"](fc)cir'kfc) 

^[l...n+l]/, X _ ^[n+l]/,x , \k) \k) \k) (AOO\ 



Note also that one can deduce physical properties directly from the recursion formulae. 
In particular, from unitarity and hermiticity relations of C^^^{k) and C[^'""'+^](/c), it is easy 
to check that these relations are also obeyed by C^^'""'~^^^{k): 

^[l...n+l] ^[l...n+l] _ J ^[l...n+l] (^^f _ ^[l...n+l] (_^^^t _ ^4^28) 

In the same way, one proves recursively that the transmission coefficients cjg" "^^^ (A;) and 
C2i"^~^^\k) depend only on the distances between the different defects, while the reflection 
coefficients obey 

^[i...n+i] _ ^[i...n+i] ^-iky, cfe'"^'^ i^) = (A:) e'^^^-^^ , 

where the 'bare' coefficients (/j;) depend only on the distances between the different 

defects. 



5 Quantum fields with defects 

This section is devoted to some aspects of quantum field theory with defects on the line. 
In order to keep the discussion as simple as possible we consider one scalar field and 
n = 2 defects. The results however can be easily extended to more general cases using 
the results of section HI 

5.1 Correlation functions 

In the Fock representation J-" {^2) of 2I2 the basic correlation functions are [17] 

{a,ik)a*ip)) = 2n{[6,,+%^{k)]6ik-p)+n,,ik)6ik + p)} , (a*(A:)a,(p)) = , (5.1) 

where i,j = 1,3 and TZ and T are given by fl3.15tl3.20|) . together with 

Tnik) = %3ik)=0, 7^13(A;) =7^3l(A;) = 0. (5.2) 

Because of (13.91 13.10|) . the correlators involving a2{k) and a2{k) can be expressed in terms 
of (15. ip . Now, using ( 12. lip one has 



{(p{tl,Xi)(p{t2,X2)) = {(Pi{tuXi)(pj{t2,X2)) = {ipj{t2,X2)(Pi{ti,Xi)) , (5.3) 

xiGDi,X2&Dj 

which can be computed for any i,j = 1, 2, 3 and have the following general form 

^— ^e-'-^'^)*- { [% + %,{k)] e'^^- + n.,{k) e'^^- } , (5.4) 
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where ti2 = ti — ^2, 2^12 = X1—X2 and Xu = Xi + X2. We recall that the dispersion relation 
is 

u{k) = Vk^ + m2 . (5.5) 

For the correlators involving z or j = 2, we use the expressions (13. 9[ 13.101) and after 
some algebra we find 

1 - 4^(-fc)i?S(A;) e^^y^- 1 - R^k)R^^l{-k) e-^^v^^ ' 
1 - Rili-k)R^iik) e^i^-s/- 1 - Riiik)R^k-k) e'^i^^- 



^ TfiJ(fc)4f(fc)e-2'^^^ 
1 - (-fc)i?i (fc) e2i%- ' 1 - (- A:)i?| (A:) e^i^^- 



T^iskfe) ^ 7j|(-fc)gj(fc)e-2i^.. 

1 - 41 (-A:)4^j (A;) e2i^^- ' 1 - 4j (fc)i?i i-k) e-2i%i2 



where we have used 'bare' coefficients to make the yj dependence explicit. 
The remaining coefficients are computed thanks to (15.31) . that is to say 

n,,{ky = TZjii-k) and %j{ky = r^i{k) , z,j = 1,2,3. 

5.2 Casimir force 

The energy density of the field ip(t, x) reads 

0{t,x) = I \{dt^r - I [^{dl^) + {dlip)^] + m'A {t,x) . (5.10) 



One can compute the expectation value {6{t,x,i)) via point-splitting from the two-point 
function (15.41) . Like in the standard Casimir effect, this expectation value diverges at 
coinciding points. This fact should not be surprising because it is the force, which is 
actually the physical observable in this context. It is given by 

d f~^°° dk 

T,{x) = —{e{t,x)) =i —kuo{k)nu{k)e^'^\ ^ = 1,2,3, (5.11) 
ax x&Di Ztt 

which determines a well-defined distribution in x. In general, Ti{x) may be singular at 
X = Ui. We will illustrate this feature below on the example of two delta defects. 
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5.3 Finite temperature 

In order to investigate the Casimir force at finite (inverse) temperature (3, we introduce 
the so called Gibbs representation Qj3{^2) of 2I2. In contrast to the Fock vacuum, the 
fundamental state of ^/3(2l2) is not annihilated by ai{k) and represents an appropriate 
idealization of a thermal bath which keeps the system in equilibrium. In Qp{^2) one 
has [17] 

{a,{k)a]{p))p = ^ _ ^^^^^^^ 27t{[5., + %,{k)] 5{k -p)+ n,,{k)5{k + p)} , (5.12) 

p-/3^(fc) 

(a*(A;)a,(p))^ = ^ _ 27r {[5,, + 7^,(fc)] 5{k - p) + n,,{-k)5{k + p)} , (5.13) 

where i,j = 1,3. It is easy to show at this point that the thermal correlation functions 
{(Pi(ti, Xi)ipj(t2, X2)) /3 are obtained by substitution 

„-/3aj(fc)+iaj(fe)ti2 i „-iw(fc)ii2 

^ e +e 

in the integrand of the respective zero-temperature expressions (15. 4p . In this way one now 
gets for the force 

= ^ 1 ^e-Sw ^^^(fe)e^''" ^ z = l,2,3, (5.15) 
showing an explicit temperature dependence. 

5.4 Charge density 

In this section we consider a complex scalar field 

ip{t,x) = ^[^^'\t,x) + iip'^'\t,x)] , (5.16) 

where both Lp^^^ and tp^"^^ are two copies of the scalar field with the same two defects 
constructed in section 3 above. The theory admits now a U{1) symmetry leading to the 
conservation of the current 

j^(t,x) = {d^ip*)ip : (t,x)- : ip*{d^ip) : (t,x)] , (5.17) 

where : ■ ■ ■ : stands for the normal product in the RT algebra 2I2. We are interested in 
the charge density 

Q{x,P) = {jo{t,x))(, (5.18) 
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of the system in the Gibbs state, defined in the previous subsection. Using the thermal 
two-point function of (/?, one obtains 



dk 



-f3uj{k) 



71 1 



-/3w(fc) 



(l+7^,,(A;)e2^'=^) 



(5.19) 



The density (15.191) vanishes in the hmit of zero temperature (/? oo): the first term in 
the integrand gives the density at finite temperature, whereas the second one generates 
the correction due to the defects. 



6 Example: two delta impurities 

We consider the case of two defects defined by the boundary conditions 

V5(i, +ya) = Vit, -ya) , 

a = 1,2, yt e 

The local refiection and transmission coefficients have the form 



(6.1) 



Rr,\k) = R^^{-k) 



k + ZyUc 



and Tl'^^{k)=n't^{-k) 



k + i^a ' 



a = 1,2. 



(6.2) 

To avoid the presence of bound states we will assume that the product /Ui /i2 is strictly 
positive. In the computation of the Casimir force we will also suppose for simplicity that 
m = 0. 



6.1 Force on a test particle between the defects 



Using eq. (15. lip , one can compute the force on a test particle when x G D2, that is to 
say yi < X < y2. We get 



TT J'2(x) 



+00 



dk- 



k' 



+ G{k) 



k — III cos(2A;(x — yi)) + /i2 cos(2A;(x — ?/2)) 
-f/x^ sin(2/c(x — yi)) + /ig sin(2A;(2; — 7/2)) 



(/ii + /i2)^ - 4/iiyU2 sm^{kyi2) + 2k fii 1^2(1^1 + ^2) sin(2A;?/i2) 



+4/ii/i2 sin {kyu] 



One can check that the integrand has no real poles, so that the divergent parts are obtain 
in the limit k 00. They can be computed analytically and we obtain: 



7lJ^2{x) 



+ 



4(a;-?/i)2 A{x-y2y 2{x - yi) 

12 



C2{x,yi,y2) - 32{x,yi,y2) , (6.3) 



where we have introduced 
C2{x,yi,y2) = 

S2ix,yi,y2) = 



dk 



+ G{k) 
G{k) 
k^ + G{k) 



- Ill cos[2A;(x - yi)] + fi2 cos[2A;(x - 1/2)] | , 
nl sin[2A;(x — yi)] + fx^ sm[2k{x — 2/2) 



The integral C2{x,yi,y2) has logarithmic singularities at x = yi, y2, x = 2yi — y2 and 
X = 2y2 — yi, which can be extracted by means of the Meijer's function 



M{x) = Gil 





OOi 



°° ,k cos(2A;x) 
dk — 



A;2 



We get 



C2{x,yi,y2) 



-fil M{x - yi) + III M{x - 2/2) - III /i2 M{x + y2 - 2yi] 
-/ii filM{x + yi - 2^/2) +'c2{x,yi,y2) 



(6.4) 



The integrals C2{x, yi, 7/2) and 32(2;, yi, 7/2) are absolutely convergent and can be computed 
numerically. We represent the force for different values of /xi and fi2 in figures [H and El 

F2(x) 
30: 

20: 

10: 



-0.7^0.^0.25 ; 0.25 

-10: 





X 



Figure 1: Force as a function of the position of the test particle x E]yi,y2[ for two delta 
impurities, with y2 = —yi = 1 and /ii = 1. 

We find a point between the defects where the force changes sign. In this point the 
test particle is in stable equilibrium. 
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-0.7^0.^0.25 

-10 

A^l=l I _20 
-30 



1=0.5 



0.250. 




X 



Figure 2: Force as a function of the position of the test particle x E]yi, for two delta 
impurities, with 1/2 = —yi = 1 and /i2 = 1- 

6.2 Force on a test particle on the left of the defects 

We perform the same kind of calculation when x E Di, that is to say x < yi. Extracting 
the poles, we get 



TT jFi (x) 



2 ul + 2fiiH2 



A;ai(A;) + ao(A;) 



/■-t-oo 

-Ciix,yi,y2) - Si{x,yi,y2) + dk 

Jo 



k^ + ^(A:) ' 



where we have introduced: 

Cl(x,2/i,?/2) = 



00 kG(k) 

dk — , <! /ii cos[2/c(a; — yi)] + /U2 cos[2/c(x — ^2)] 







Si(x,yi,?/2) 



k^ + G{k) 
G{k) 



dk 

k' + G{k) 



/ii sin[2/c(a; - yi)] + (//2 + 2/xi/i2) sin[2A;(x - 2/2)]} , 



ai{k) = 2/ii/i2|/i2 cos[2A;(x - - /i2 cos[2fc(x - 2/2)] + /^i sin[2fc(x - 2/1)] sin(2A;yi2)| 
ao(A;) = 2/i^/i2 sin[2A;(x — yi)] [l — cos(2A; 1/12)] . 
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Again, one can extract from Ci(x, 2/2) the Meijer's contribution 

Ci{x,yi,y2) = -fj.lM{x -yi) + fi2{if^i + f^2Y + 2fil) M{x -y2) 
-Hi nl M{x + yi- 2y2) +'ci(x, yi, 2/2) • 



(6.5) 



The remaining integrals are absolutely convergent and can be computed numerically. We 
represent the force in figures [3] and HI which show that it is repulsive close to the defect 
at yi, but for x < yi there are also domains of attraction. 



F1(X) 

2r 





Figure 3: Force as a function of the position of the test particle x < yi for two delta 
defects, with 2/2 = —yi = 1 and /ii = 1. 




^l=0 . 7 
^1=0. 5 
ij1=0.3. 



Figure 4: Force as a function of the position of the test particle x < yi for two delta 
defects, with y2 = —yi = 1 and /i2 = 1. 
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6.3 Force at the defect 



As already mentioned, the force is a distribution with singularities at the positions x = 
yi, y2 of the defects. To compute the force acting for instance at the defect in yi, one 
needs to sum the two contributions computed in sections 16.11 and 16.21 and take the limit 
X ^ yi, getting rid of the divergences. This is achieved by taking the symmetric limit: 

^12 = hm - 5) + T^iyi + . 

One obtains the force as a function of the distance between the two defects: 

1^2 , /i2(/i2 + /il) , , , , G{k) 



12 



2Z/?2 

-2/ii fij 



yi2 

kdk 



i) f Gik) 

- - 2/i2(/ii + 112) dk ^ ^^^^ sin(2/c 2/12) 



,0 k^ + G{k) 

r+oo ^3 



k^ + G{k) 

(yUi + /i2) k sin(4?/i2) + 4/ii^2 cos(2/c?/i2) sin^(/c?/i2) 
yUi/i2 cos(4/c?/i2) + + 1^1 + ^1^2) cos(2fc?/i2) 



Using the Meijer's function, one finds for the last integral 

-2/ii/z^ M{2yi2) - 2^2 {fJ-l + fJ'l + /^i/^2) M{yi2) + Ci2(?/i2, fJ-i, ^^2) , 

where £12(^/12, /^i, /^2) is absolutely convergent. Computing the integrals numerically, one 
gets the figures O and [Hi 



F12 




yl2 



^2=l 
^j2=0.5 
2=0. 35 



Figure 5: Force at the defect for two delta defects, as a function of the distance between 
them for /ii = 1. 
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F12 




^1=1 

^1=0.5 

^.1=0. 35 



Figure 6: Force at the defect for two delta defects, as a function of the distance between 
them for /i2 = 1- 



6.4 Charge density 



The charged density for a complex field and two delta-defects is obtained from fl5.19p . 
The integrals can be computed numerically, the mass m and the temperature ^ playing 
the role of regulators at k = and k = oo respectively. We get 



pjix,P) = pjix,(3) 



dk 



-I3u){k) 



dk 



-f3uj{k) 



Xj[X) 



1 - e-z^-'W k^ + G{k) 



J = 1,2, 



with 



Mk) 



XiiX 



2 sm[2k{x — Hi)] — kp2 sin[2A;(x — j/2)]| k^ 
+2 cos\2k{x — Ui)] + p\ cos[2A;(x — 2/2)] | k^ , 
2 sm[2k{x — 2/1)] + /U2 sin[2A;(a; — ?/2)] | k^ 
-2 cos[2fc(x - yi)] + /U2(2/ii + ^2) cos[2/i;(x - 2/2)] j^:^ 
-4 yu,i/i2 sm{2kyi2) \ fJ-i cos[2/i;(x - yi)] + ^2 cos[fc(2a; - yi - 2/2)] \k 



-Afijpl 



cos{2kyi2) cos[2k{x — yi)] 



and a similar expression for r3(A;). We remind that uj{k) = y/k^ + m?. 

When m = 0, there is a divergence at A; = reflecting the possibility of Bose condensa- 
tion. When m 7^ (and /3 7^ 0), the integral defining pj{x, P), j = 1 or 3, converges. The 
same is true for p2{x, j3) and generic values of the parameters /xi, //2, yi and ?/2- However, 
when the relation 

Pi + P2 = 2/U1/U2 {y2 - yi) , 
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is fulfilled, P2{x,(3) diverges (see figure [7]). This relation corresponds to a change in the 
behavior of G{k) around k = 0: 



G{k) (Afi + /i2 - 2^1 /i2 2/21)^ + o(A;'^) , 

k = 

and could be interpreted as a kind of resonance effect. One can check that, at this point 
in the parameter space, the integrals entering in pj{x,f3), j = 1 or 3, as well as in the 
force studied in the previous sections, still converge. 
The plots of pj{x,(3) are presented in figures IHlfTOl 



-5000 
-10000 
-15000 
-20000 
-25000 



0.2 



.4 0.6 



x=0 

x=-0.5 



Figure 7: Density between two delta defects, as a function of ^1 for /i2 = 1, 1/2 = 1 = —yi 
and P = 1. 



P~l(X) 




Figure 8: Density on the left of two delta defects, as a function of x for /^i = 1, j/2 = 1 = 
—yi and [3 = 1 (large distance behavior). 
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/ii = 0.35 111 = 0.2 



Figure 10: Density between two delta defects, as a function of x for /i2 = 1 and (3=1. 
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7 Conclusion 



In the present paper we generalized the RT algebra framework [14]- [16] for dealing with 
one defect to any number of defects. More precisely, we considered a scalar field, which 
interacts with multiple defects on the line and freely propagates away of them. We studied 
the most general point-like defects, parametrizing the interaction by local reflection and 
transmission coefficients. As expected, the algebraic formulation considerably simplifies 
the boundary value problem at hand and is very efficient for deriving the correlation 
functions both at zero and finite temperature. Applied to the Casimir effect, the algebraic 
technique leads to a nice and compact expression (in terms of the reflection coefficients) 
for the force. As an explicit illustration, we derived the Casimir force for two delta-type 
defects. It turns out that both the intensity and the direction of the force depend on 
the position. We established also the charge density distribution of a complex scalar fleld 
interacting with the defects. 

Our framework can be generalized in several directions. One can apply it to the case 
where the fleld ip has > 1 internal degrees of freedom. Then, each of the operators 
ai(/c), z = 1, . . . ,n, becomes a vector of length A^, while the 'coefficients' C|°', i,j = 1,2 
are now Nx N matrices. The above results remain valid if these matrices commute among 
them. 

Another very attractive application concerns quantum graphs [21,22], which are net- 
works of wires connected at nodes. Each node is characterized by a scattering matrix [23] 
and can be viewed as a defect [24]. The RT algebra approach has been already applied [25] 
to the study of quantum fields on graphs with one vertex (star graphs). Using the results 
of this paper, one can formulate quantum field theory on a generic quantum graph and 
address in this framework the interesting problems of conductance and vacuum energy. 
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